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PACS numbers; 87.10. +e; 87.16.-b 

In their elegant work^ Aurell and Sneppen developed 
a framework to quantify the robustness of a biological 
state. They have successfully applied their formulation 
to the robustness of the epigenetic state of the wild type 
phage A, but concluded, and supported by their fur- 
ther numerical study^, that their model is not robust 
against various mutations in the gene regulatory network 
of phage A. They called this disagreement between their 
mathematical modeling and biological experiments the 
stability puzzle. We agree with their assessment that the 
robustness can be viewed as a first exit problem and do 
not dispute the possibility that certain unknown biology 
may be responsible for their discrepancy. In this com- 
ment we would like, however, to call the attention to two 
important aspects, one on biology and one on modeling, 
overlooked by Aurell and Sneppen, and to argue that 
the including of those two aspects may be sufficient to 
understand the stability puzzle. 

Biological aspect: in vivo vs. in vitro. It is well 
known that so far most molecular parameters in biology 
are measured in vitro, that is, not in the native chem- 
ical and biological environments of the living organism 
{in vivo). This is same for phage A, one of most well 
studied organisms.^ Tremendous efforts have been made 
to ensure the in vitro conditions as close as those of in 
vivo. Subsequently, it has been generally concluded that 
though there is indeed an in vivo and in vitro difference, 
this difference is small and normally not exceeds 30% of 
the measured in vitro value. ^ Such a small difference, for 
the case of affinities in DNA-protein binding dominated 
the dynamics of gene regulatory networks, is about a few 
kcal/mol. Regarding to the usual qualitative nature of 
biological data where measurements for the same quan- 
tity are common to be differed by a factor of order unity, 
one would be satisfied to have such in vitro values. Un- 
fortunately, for the stability as a first exit problem, the 
escape rate can be exponentially sensitive to some pa- 
rameters. Hence 10% change in a critical parameter can 
result in orders of magnitudes change in the escape rate. 
In physics and chemistry such examples are too abun- 
dant to enumerate here. This comes our first observation 
to understand the stability of the biological states: One 
must carefully take the in vivo and in vitro difference into 
account for dominant molecular parameters in the mod- 
eling. This has been attempted but not taken seriously 
in Ref.[l] and [2]. 

Modeling aspect: the nature of noise. As it is 
clear from the mathematical modeling, the noise strength 



equivalent to temperature plays the most important role, 
as well exemplified by the classical Arrhenius law. We 
found two important features on noise which have not be 
taken cared of properly in Ref. [1] , the diffusion approxi- 
mation for the intrinsic noise and the presence of extrinsic 
noise. First, Aurell and Sneppen modeled the noise as it 
directly comes from the chemical reaction rates. Indeed, 
the chemical reactions are noisy processes. What used 
by Aurell and Sneppen is the standard diffusion approx- 
imation to obtain the noise strength, called the intrinsic 
noise. ^ This approximation is good near a stable fixed 
point, but can be very problematic when applied to rela- 
tive stability problem. Orders of magnitudes in error may 
be embedded in this diffusion approximation when cal- 
culating the escape rate.^ More precise approximation is 
needed to get a better description of the intrinsic noise,'' 
which we believe can be done with an improved diffu- 
sion approximation. The second feature is more seri- 
ously, regarding to the extrinsic noise: There are other 
noisy processes beyond those dictated by reaction rates 
of the deterministic part of the modeling. In fact, bio- 
logical experimental data already indicated the existence 
of the extrinsic noise. ^ The extrinsic noise is not in the 
modeling of Aurell and Sneppen, which may well account 
for their inability to solve the stability puzzle. 

With the consideration of above two biological and 
modeling aspects, the in vivo and in vitro difference and 
the proper modeling of the noise, can one understand the 
stability puzzle in phage A? It is evident that those two 
aspects will give enough difference in the numerical out- 
come of the modeling. One should not be surprised to 
find the positive answer. 
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